Numerical Study of Photo-Induced Dynamics in Double-Exchange Model 
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Photo-induced spin and charge dynamics in double-exchange model are numerically studied. 
The Lanczos method and the density-matrix renormalization-group method are applied to one- 
dimensional finite-size clusters. By photon irradiation in a charge ordered (CO) insulator associated 
with antiferromagnetic (AFM) correlation, both the CO and AFM correlations collapse rapidly, and 
appearances of new peaks inside of an insulating gap are observed in the optical spectra and the 
one-particle excitation spectra. Time evolutions of the spin correlation and the in-gap state are 
correlated with each other, and are governed by the transfer integral of conduction electrons. Re- 
sults are interpreted by the charge kink/anti-kink picture and their effective motions which depend 
on the localized spin correlation. Pump-photon density dependence of spin and charge dynamics 
are also studied. Roles of spin degree of freedom are remarkable in a case of weak photon density. 
Implications of the numerical results for the pump-probe experiments in perovskite manganites are 
discussed. 

PACS numbers: 71.30.-l~h, 78.47.J-, 78.20.Bh, 71.10.-w 
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I. INTRODUCTION 



Coherent controlling of electronic and structural 
phases by shining light has been one of the attrac- 
tive themes not only in fundamental condensed-matter 
physics but also in electric and communication technolo- 
gies. Drastic state-change by photon irradiation, often 
termed photo-induced phase transition, is ubiquitously 
seen in a variety of materials In particular, a num- 
ber of studies in photo-irradiation effects have been done 
in correlated electron systems, such as transition-metal 
oxides, low-dimensional organic salts and others. This 
is because i) cooperative and ultra-fast changes of states 
by light occur due to strong electron correlation, and ii) 
coupling between multiple degrees of freedom of electron, 
i.e. spin, charge and orbital, exhibits a variety of photo- 
induced phenomena. In recent rapid progresses of exper- 
imental techniques, transient and non-equilibrium states 
induced by light irradiation have been examined by sev- 
eral kinds of time-resolved scattering and spectroscopy 
experiments, for example photoemission spectroscopy^^ 
x-ray diffraction;^!^ electron diffraction^ and so on in 
addition to the conventional optical pump-probe spec- 
troscopy. 

Perovskite manganites Ri-xAxMnO^ {R: a rare-earth 
ion, A: an alkaline-earth ion) is one of the well studied 
correlated electron materials in a viewpoint of ultra-fast 
photo-induced phenomena. A nominal valence of a Mn 
ion is 3 4- x and the electron configuration of the 3d or- 
bitals is One of the doubly degenerate Cg orbitals 

is occupied by 1 — a; electron and the t2g orbitals are 
occupied by three electrons with parallel spin. One re- 
markable property in manganite is strong phase compe- 
tition. Around x = 0.5, a subtle energy balance controls 
stability of two electronic phases; charge ordered (CO) 
insulator associated with antiferromagnetic (AFM) or- 



der and ferromagnetic (FM) metal<i Colossal magnetore- 
sistance (CMR) effect is one of the examples of drastic 
phase change between the two phases by applying exter- 
nal ficldj^ Electronic and structural phase controlling is 
also performed by light i^"— By photo- irradiation in CO 
insulating phase, the optical absorption spectra around 
the charge-transfer excitation shift to lower energy iii Af- 
ter 10ns, the spectra are almost recovered. These exper- 
imental results imply a generation of a transient metallic 
state by photo-irradiation and its fast relaxation within 
10ns. The transient magnetic property was examined by 
the magneto-optical Kerr spectroscopy. A Kerr rotation 
appears by photo irradiation and its angle gradually in- 
creases within Ipsii^ii^iii That is to say, spin and charge 
structures are changed cooperatively by pump photon 
irradiation. The excitation density dependence of the 
photo-induced state was also examined in the pump- 
probe expcrimentSFi^ The results provide a clue to re- 
solve a stability of the photo-induced FM metal and its 
relaxation mechanism. 

Theoretical study of manganites has been done mainly 
from the view point of origin of CMR. The double- 
exchange (DE) model is not only studied as a model for 
manganites, but also recognized to be a standard theo- 
retical model in solid state physics<^"— In this model, 
conduction electrons corresponding to the Cg electrons 
couple ferromagnetically with localized spins for the t2g 
spins. On the contrary, transient and non-equilibrium 
properties in the DE model have not been examined well 
so far. 

In this paper, motivated from the optical pump-probe 
experiments in perovskite maganites, we present a nu- 
merical study of the photo-induced dynamics in the one- 
dimensional DE model. Time dependences of the static 
and dynamical quantities are calculated by applying the 
Lanczos method and the density-matrix renormalization- 
group (DMRG) method in finite size clusters. By photo- 



2 



irradiation in the CO and AFM insulating phase, both 
the CO and AFM correlations rapidly collapse. New peak 
structures appear inside of the gap in the optical absorp- 
tion spectra and the one-particle excitation spectra, and 
grow up with time. It is shown that time dependences 
of the spin correlation and the new peak structures are 
scaled by a universal curve. This implies strong coupling 
between the charge and spin sectors in transient photo- 
excited state. We also study the pump-photon density 
dependence of the photo-induced dynamics. Roles of spin 
degree are remarkable in a region of weak photon den- 
sity in comparison with charge degree. Implications of 
the present numerical calculations for the photo-induced 
phenomena in manganites are discussed. 

In Sect, nil the extended DE model and a formulation 
of the photo-excited states are presented. The numerical 
methods applied to one-dimensional clusters are intro- 
duced in Sect. IIIIl The ground state properties before 
photo-irradiation are briefly introduced in Sect. IIVI The 
main part in this paper is Sect. Ivl where the numerical 
results for the photo-induced spin and charge dynamics 
are presented. The pump-photon density dependences of 
the photo-induced dynamics are shown in Sect. I VII Sec- 
tion [Vn] is devoted to discussion and concluding remarks. 
Details of an effective model and analytical results in the 
spin-less V -~t model arc given in Appendix A. A part of 
this paper was briefly presented in Ref. 



II. MODEL AND FORMULATION 

We introduce the one-dimensional extended DE model 
to examine both the CO insulating state associated with 
the AFM order and the FM metallic state. A model 
Hamiltonian is given by 

■Ho = -at ^ (^cL^jcr + H.c.^ + X] "»t"»4- 

(jj)cr ■( 

where Cia is the annihilation operator for a conduction 
electron with spin a{=^, 4) at site i, and Si is a spin op- 
erator for the localized spin. We introduce the number 
operator rii = ^i'^ ~ '^L'^io' ^^'^ ^P^^ oper- 
ator = (1/2) X^a/j cJq i^)af} "^ith the Pauh matri- 
ces <T for the conduction electrons. The first term in 
Eq. ([T|) is for the electron hopping with amplitude at be- 
tween the nearest neighboring (NN) sites. The second 
and third terms represent the on-site Coulomb interac- 
tion U, and the NN Coulomb interaction V for the con- 
duction electrons, respectively. We consider the Hund 
coupling Jh{> 0) between the conduction electron and 
the localized spin, and the AFM superexchange interac- 
tion Js{> 0) between the NN localized spins. We take 
t = 1 as a unit of energy, and a magnitude of the transfer 
integral is changed by changing the parameter a. For 



simplicity, we assume a single orbital for a conduction 
electron and 5 = 1/2 for a localized spin. 

The vector potential for the pump photon at time r is 
given as Apuinp(r) = ^pump(''')'2 where the z axis is taken 
to be parallel to the chain direction. The interaction 
between the conduction electrons and the pump photon 
is introduced as the Peirles phase: 



n'{r) 



-at 4a 



-i f Ap^„-,p{T)-dr 



H.c. 



(2) 



where a subtraction in the parenthesis corresponds to 
the first term in Eq. ([1]). In a case of weak pump- photon 
density, this Hamiltonian is reduced into the following 
form, 



(3) 



up to the order of 0(v4pu„ip). We introduce the current 
operator defined by j = iatJ2(ij)a (cLcjcr - H.c^. By 
using the first-order time-dependent perturbation theory 
with respect to 'H'(t), the wave function at time r is 
obtained as 



|0> 



Cnir) \n) , 



(4) 



where a coefficient is given by 



c„(r)=ze-^^"-(n|j|0) 



dr'A 



pump 



(r')e 



-i{Eo-E„)T' 



(5) 



and \n) is the n-th eigen-state of Hq with the energy En- 
We assume that the system is in the ground state |0) at 
T = —oo, and focus on the second term in Eq. ^ which 
represents the one-photon absorbed state. From now on, 
this is termed \4'{t)) as 



|0(r))^^c„(r)|n) 



(6) 



The vector potential for the pump-photon is assumed to 
be a damped oscillator form as 



^pump(T) = Ai exp (-7o|r| - iujor) 



(7) 



where is a frequency and 70 is a damping factor. We 
impose a condition of t ^ which implies that the 
pump photon is fully damped at time t of our interest. 
Then, the coefficient in Eq. ([5]) is obtained as 



= -2iAie-'^"^ {n\j |0) Im 



1 



^0 — En 



En 



no 



(8) 
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probe 




FIG. 1: A schematic view of the wave packets for the pump 
and probe light intensities as a function of time. 



By inserting this expression, we have the wave function 
at time t as 



|0(r)) = -2*^1^6- 



X Im 



1 



ujQ - En + En+ ijo 



\n) {n\j\0) 
1 



Wo - Ho + -Eo + 470 



J |0) 



(9) 



where (0|j|0) = is used, 
photon absorbed state by 

1 



Here we redefine the one- 



(r)) = -^e-'«°nm 



^0 ^ Ho 



En 



«7o 



J |0), 



(10) 



with a normalization factor TV. 

Photo-excited states are monitored by calculating the 
several kinds of static and dynamical quantities. The 
transient excitation spectra are formulated based on the 
liner-response theory. We explain, as an example, a for- 
mulation of the transient optical absorption spectra. The 
vector potential for the probe photon is set to be parallel 
to the z direction and has a damped oscillator form as 
^probo(T) = Ap^ohc{T)z and 



^probo(T) = A2 exp {-j\t - td \ ~ iujT) . 



(11) 



with a frequency uj, the center of an envelope tjj, and 
a damping factor 7. We assume that the pump pho- 
ton is fully damped at time when the probe photon 
comes in, and there are no interference between the 
pump and probe photons. This is given by the condition 
{td — 7^^) ^7o^^. A schematic picture for the pump- 
photon and probe-photon wave packets is shown in Fig.[T] 
The interaction Hamiltonian between the probe photon 
and the conduction electrons at time r is given by 



H"{t) = -J^probc(T), 



(12) 



within the first order of ^probo (''')• When the probe pho- 
ton is taken into account, the electronic wave function at 
time r (» 7(7^) is given by 



|^(r)) = C/(T,ro)|(/)(To)). 



(13) 



Here, \(f>{To)) is the wave function at time tq where a 
condition To ~ ^ tq '3> is satisfied and is given 
by Eq. (fTO|l . The time-evolution operator is defined by 



with 



UiT,To) = Uo{t,to)Ui{t,to), 
Uo{t, Tq) = exp [-iHo (t - tq)] 



and 

Ui{t,to) = Pexp 



i / dr' j{T' , To) Apy-ohcir') 



(14) 
(15) 

(16) 



where P is the time-ordered operator, and 

j(t',to) = Uq{t' ,TQ)jUo{T' ,to) is the interaction 

representation of the current operator. The expectation 
value of the current operator at time r is given as 

(j)(r)^(^(r)| J |^(r)) 

= [ dr' (</>(ro)| [j{T',To),j{T,To)] |0(to)) Ap,.obc(T') 

Vobc;- (17) 



J To 

0«obe)- 



Here we use a relation (0(to)| j(r, tq) \(j>{To)) = 0, because 
|(/)(t)) is an eigen state of parity. By inserting the com- 
plete set into the integrand and integrating out r' under 
the condition of (tu — tq) ^ 7~^, we obtain 

(j) (t) = X{^,T)Ap,ohc{T) + Xo{oJ,T)Aprohc{TD), (18) 



with 



x(w,t) = -^<('r)Q(r)j„„j„ 

nml 

1 



X 



OJ - Em +Ei±ij OJ - En + Em ± ij 



(19) 



and 



Xq{^, t) = -e{T - TD)2ij} jnmjml 



ninl 



<iT)ci{TD)e- 



iE„,{T-ru) 



{uj - Em + EiY +7' (w - E„ + E„ 



(20) 



We introduce c„(t) = (n|0(T)) which is proportional to 
c„(r) in Eq. ([5]), and j„m = {n\j\'m). The first and 
second terms in Eq. (|18p represent the induced current 
which is proportional to the probe-photon vector poten- 
tial Api-obc('''), and the oscillating part which remains to 
be finite at r = +cxd, respectively. The positive and neg- 
ative signs of 17 in the denominators of Eq. (jl9p are for 
the cases of r < td and t > td, respectively. We focus 
on the first term in Eq. ^TE\i with the positive sign of 
17 in x(a;,r), from now on. The first and second terms 
in x(a;,T) represents the photon-absorption and photon- 
emission processes, respectively. The optical absorption 



4 



spectra is defined by the imaginary part of the first term 
given as 

a{u:,r) = _ ±i^y ^ilhi^lll^ ^ (21) 
^ ' Lit ^w-E^ + Ei+ij ^ ^ 

nml 

where L is the system size. 

In the similar way, the one-particle excitation spectra 
are obtained as a sum of the electron and hole parts: 
A{q, uj) = A^°^''{q, uj) + uj) with 



and 



^ <(r)c;(r) (?7.| Cq„ \m) (m| \l) 



nralcr 



Im- 



-7r(iLj - Em + Ei+ ij) 



(22) 



nralcT 



Ini 



<{t)ci{t) {n\ c\„ \m) (m| C-q„ \l) 



-tt{-uj - Ejn + Ei+ 17) 



(23) 



We also obtain the imaginary part of the dynamical 
charge susceptibility as 

= ~-^cI{t)ci{t) {n\ng \m) {m\n_q \l) 



X Im 



nml 

1 



cj — Em + El + i'-f UJ 



(24) 



and that of the dynamical spin susceptibility for the lo- 
calized spin as 



S{q,L 



X Im 



y!<Mcz(T) {n\ Sg \m) {m\ S^g \l) 



nml 

1 



oj — Em + El + i-f CJ — E„ 



E„ 



■ 27 



(25) 



where Cg^, ng, and Sg are the Fourier transforms of Cia, 
rii, and Si, respectively. 

In the numerical calculation, to reduce the computer 
resource, En and Ei in the denominator of Eq. (|19p are 
replaced by the expectation value for T-Lq with respect to 
the wave function after the pump irradiation, i.e. Ep = 
{4>{t) \ Tio |(/)(r)). Then the response function is rewritten 
as 



X(w,r) ~ - (0(r)| j 
1 



UJ - Ho + Ep + i'y UJ - Ep + Ho + 



(26) 



The same approximation is applied to the optical absorp- 
tion spectra and the one-particle excitation spectra. This 
approximation is reasonable in the present case where the 
pump-photon energy is tuned around a narrow energy 
region of wq ± 7o with 70 <Si ujq. We have numerically 
confirmed in the calculation of a(uj) that this approx- 
imation reproduces quantitatively the results based on 



the exact expression in Eq. ((2T|) . In order to unify the 
expressions before and after pumping, we introduce the 
following expression for the optical absorption spectra: 



a(uj) = Im / 7 ; ^ 

^ ' Ln Y UJ -no + (Ho) +h 



(27) 



A bracket (• ■ • ) implies the expectation value with re- 
spect to |0) before pumping, and that with respect to 
|(/)(t)) after pumping. 

III. METHOD 

In order to analyze the photo-excited states in the gen- 
eralized DE model, the ED method based on the Lanczos 
algorithm and the DMRG method^^"— are applied into 
one-dimensional finite-size clusters. 

In the ED method, one-photon absorbed state at r = 
is calculated by inserting a subset {\n)} with finite dimen- 
sion Ml, instead of the complete set in Eq. pH)) 
as 



10(7 



Ml 

0))^-^;^Im 



J|0) 



i^O — En + Eq 



no 



(28) 



This set of the wave functions is obtained by the Lanczos 
procedure with the Mi steps from a trial function j |0), 
and Efi's are the corresponding eigen energies. We take 
Ml = 300—400 which is enough for the optical absorption 
spectra before pumping. This reduction technique is also 
used to calculate the time evolution of the wave function. 
The wave function at time t+St is obtained from |0(r))^ 



\Ht + St)) 



Ah 



n) imr)) , (29) 



where a subset with a dimension is obtained by the 
Lanczos procedure for the Hamiltonian "Hq from the trial 
function |0(r)). Wc take A/2 = 20 and 5t = IQ-'^/t - 
/t which are enough for the several time-dependent 
quantities. We checked that the total energy remains to 
be a constant within a numerical error. 

In the DMRG method, wc use the muhi-target DMRG 
algorithm to obtain the excitation spectra and the time 
evolution. The density matrix for the system block is set 
to be 



aj 



(30) 



where p*^"' is the weight for the a-th target with a rela- 
tion X^Q^''""'' ~ 1' j represent the states in the 
system and environment blocks, respectively, and I/;'-"-' 
is the wave function of the a-th target state. For the 
excitation spectra with respect to the operator O before 
pumping, we take the following target states: 



|0),O|0),. 



1 



Ho + Eq + 



O\0) 



(31) 
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For example, we take O = j for the optical absorption 
spectra. In the calculation of the correction vector in 
Eq. (|3T|) . we use the reduced basis set obtained by the 
Lanczos procedure from the trial wave function O |0). To 
obtain the wave function \4>{t)) at time r after pumping, 
we take the following target states: 



{|O),j|O),|0(r))}, 



where we define 



mr)) 



^ Ms 



^^Im 



(32) 

\n). 
(33) 



Here, {|7"i)} is the reduced basis set with dimension 
M3 = 300 — 400 obtained by the Lanczos procedure in 
"Ho from the trial function j |0), and iS^'s arc the corre- 
sponding eigen energies. The total energy remains to be 
a constant within numerical errors. In the calculation of 
the transient excitation spectra at time t, the following 
target states are adopted, 



|O),j|O),|0(r)),O|0(r)), 



1 



Er, 



17 



-o\m) 



(34) 



The truncation number is chosen to be to = 300 in most 
of the calculations, and to = 800 in some cases. 

To remove boundary effects, system size L is set to 
be odd, and a number of the conduction electrons are 
N ~ (L + l)/2 which correspond to the quarter filling 
in the open-boundary condition. We take L = 9 and 13, 
and confirmed no qualitative differences between the two 
results. We mainly show the results in L = 9 and iV = 5. 



IV. INITIAL STATE 

First we show the electronic structure before photo- 
irradiation. In Fig. [2ja), the phase diagram is presented 
in the plane of V and J5. This is obtained by the charge 
correlation function defined by 



2 ^ 

N{q) = ^ ^ _^ sin qri singr^ {SuiSrij) , 



(35) 
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FIG. 2: (color online) (a) Ground-state phase diagram. Dots 
represent the calculated points. Abbreviations CO, M, FM, 
AFM and IL imply the charge ordered phase, the metal- 
lic phase, the ferromagnetic phase, the antiferromagnetic 
phase and the island phase, respectively. The inset shows a 
schematic picture for the spin and charge configurations in the 
CO-AFM phase. A circle represents the parameters adopted 
in the numerical calculation in Sects. |V] and IVII (b) Charge 
correlation function N{q), and (c) spin correlation function 
S{q) for the localized spins. Symbols I-V correspond to the 
parameters in (a). Numerical data of S{q) in I and II are 
almost overlapped. Other parameters are chosen to be L = 9, 
U = 20t and Jh = 15t. 



with the charge fluctuation operator dui = m — N/L, 
and the spin correlation function for the localized spin 
defined by 



Siq) 



— ^ sin qri sin qr^ {Si ■ Sj) . (36) 



The correlation functions in each phase are shown in 
Fig. mb) and Fig. HJc). Because of the open-boundary 
condition adopted in the present cluster, we introduce the 



quasi-momenturo^L in a potential well of width L defined 

by 

nn 



q 



L + r 



(37) 



with an integer n{= 1, 2, • ■ • , L), and the momentum rep- 
resentation of a local operator Oi as 



sm qri 



(38) 



6 



0.15 




FIG. 3: Optical absorption spectra q(cij). Parameters are 
chosen to be 7 = 0.2t and L = 9. 



which takes zero at i = and L + 1. It is shown in the 
figures that, in the CO phase, N{q) takes a sharp peak 
at g = TT. There are three CO phases; 1) the CO-FM 
phase, 2) the CO island (CO-IL) phase where S{q) has a 
weak peak, and 3) the CO-AFM phase. In the CO-AFM 
phase of our present interest, total spin-quantum number 
defined by = 0t + is (i + 3) /4, and a correla- 
tion function between a localized spin and a conduction 
electron at the same site is N~^J2i{'^i ' = 0.247. 
These imply a ferrimagnetic configuration as shown in 
the inset of Fig. [2ja) where the conduction electrons and 
localized spins form the spin-triplet states at every other 
sites. The obtained phase diagram is qualitatively con- 
sistent with the previous result in Ref. |32| 

The photo-induced spin and charge dynamics are ex- 
amined in the CO-AFM phase near the phase boundary. 
A set of the parameter values are U = 20t, V = 5t, Jh = 
15t, Js = 0.3< which are marked by a circle in Fig. [IJa). 
These values are somewhat larger than the realistic values 
for the pcrovskite manganite, but are required to realize 
the CO-AFM phase in one-dimensional clusters. 

The optical absorption spectra calculated in this pa- 
rameter set is shown in Fig. [3l In naive sense, these 
optical spectra correspond to the excitations from the al- 
ternate charge configuration • • ■ 1010101 • • • in a chain to 
the configuration • ■ • 1011001 • • ■ where 1 and represent 
the electron occupied and unoccupied sites, respectively. 
In the next section, we focus on the photo-induced dy- 
namics where the pump-photon energy is tuned around 
the edge of these spectra. 



PHOTO INDUCED DYNAMICS 




FIG. 4: (color online) (a) Charge correlation functions, and 
(b) spin correlation functions of localized spins at various 
times, (c) Time dependence of the difference energy for each 
term in the Hamiltonian AEx- The NN Coulomb-interaction 
term is scaled by 1/10. 



at the lowest peak in aiuj) in Fig. |31 i.e. = 3.8t with 
a damping constant 70 — QAt. 



A. Static Correlations 

The charge and spin correlation functions for several 
times are shown in Figs.|4l^a) and|4l^b), respectively. Just 
after the photo- irradiation, sharp peaks in N{q) and S{q) 
around q = ir are suppressed. After rt = 3, N(q) is 
almost independent of time, but S{q) in large (small) q 
regions decreases (increases) weakly with time. We also 
show the time dependence of the expectation value for 
each term in the Hamiltonian in Eq. ([T]). We define 



In this section, we introduce the numerical results for 
the photo-induced dynamics. The pump energy is tuned 



-at 



E 



4aCj<J + H.C. 



(39) 
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Eh = -Jff^(5,-s,;), 



E. 



's - Js^{Si ■ Sj), 



(40) 



(41) 



(42) 



(43) 



and differences between the expectation values at time 
t(> 0) and those before photo-irradiation defined by 



AEx{T)=ExiT) 



Ex : 



(44) 



for X = it, U, V, H, S). Brackets (• • • ) in Eqs. (I2S1)-(I1S1) 
represent the expectations in terms of \(I>{t)) and |0) for 
Ex{t) and E^^, respectively. Results are plotted in 
Fig. |4] (c). Just after the photo- irradiation, a large in- 
creasing of the Coulomb interaction term AEy implies a 
melting of CO. This term does not show remarkable time 
dependence furthermore, and the initial CO state is not 
restored within this time scale. Remarkable time depen- 
dence is also seen in the kinetic energy term AEt , and is 
almost compensated by that in AEs- That is, the energy 
flows from the conduction electrons to the localized spins. 
Neither remarkable change by the photo-irradiation nor 
large time dependence are seen in AEh- The local spin- 
triplet states between the conduction electrons and the 
localized spins are maintained, because the pump-photon 
energy is much smaller than the Hund coupling. 



B. Optical Absorption Spectra 




FIG. 5: (color online) Optical absorption spectra at various 
times. Broken line is for the dynamical correlation function 
for the stress-tensor operator before pumping. System size is 
chosen to be L = 9. 
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The optical absorption spectra before and after pump- 
ing are presented in Fig. [SJ We also plot the dynamical 
correlation function for the stress-tensor operator before 
pumping. This is defined by 

= - -^Im (0| 6tj n, ^- |0) , (45) 

where Sr] is the stress-tensor operator given by Srj ~ ij — 
(0| 77 10) and 



H.c. 



(46) 



A parity in the final states in /3(cj) is the same with a par- 
ity in both the ground state |0) and the final states in the 
transient-optical spectra a(w) after photo-irradiation, 
but is different from a parity in the one-photon absorbed 
state \4>{t)) in Eq. (|10p . After the photo-irradiation, a 
new peak appears inside of the optical gap and grows up 



FIG. 6: (color online) Optical absorption spectra in the spin- 
less V — t model. Broken line is for the dynamical correla- 
tion function for the stress-tensor operator before pumping. 
Numerical data for rt =3, 6, 10, 15 are almost overlapped. 
Parameters are chosen to be V/t = 5, L = 9 and N = 5. 



with time. A position of this peak around w ^ 1.5t cor- 
responds to the energy difference between the two peaks 
denoted by A in a{uj) and B in (3(uj) before pumping 
shown in Fig. [S) This in-gap spectral weight is attributed 
to the excitations between the two, and corresponds to 
the Drude component in the thermodynamic limit as ex- 
plained later. It is checked that a small weight around 
uj = 0.2t vanishes in the large limit of V/t. 

In order to understand, the photo-induced in-gap state 
in a{uj) in more detail, we focus on the charge degree 
of freedom and study the spin-less V ~ t model in one- 
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dimensional chain defined by 



vt 



H.c. 



(ij) 



(47) 



Here di is the annihilation operator of a spin-less fermion 
at site i and rii = d\di is the number operator. We 
consider the half-filling case where the average number of 
fermion per site is 0.5. The Lanczos algorithm is applied 
into the one-dimensional cluster with the open-boundary 
condition. The optical absorption spectra are shown in 
Fig. [6] By photo-irradiation where the pump photon 
energy is tuned at = 3.2f with a dumping factor 70 = 
0.4i, a new peak appears inside of the optical gap. This 
component is attributed to the excitation between the 
peak A in a(a;) and the peak B in /3(a;) in Fig. [51 as well 
as the in-gap component in the DE model shown in Fig.[S] 
Therefore, the in-gap component shown in the DE model 
is attributed mainly to the charge degree of freedom and 
its mechanism is able to be examined within the V — t 
model. On the other hand, the time dependence of the in- 
gap band in the DE model is qualitatively different from 
that in the V — t model where the new peak is almost 
time-independent after photo-irradiation. Spin degree of 
freedom in the DE model plays main roles on the time- 
dependence of the in-gap component. 

A charge configuration before pumping and that in the 
charge-excited state in the V ~t model are schematically 
represented in a one-dimensional chain as •••lOlOlOl--- 
and • • • 1011001 • • • , respectively. Here, 1 and imply the 
electron occupied and unoccupied sites, respectively. In 
the photo-excited states, there are one NN electron pair 
represented by 11 and one hole pair 00 which are equiv- 
alent to a kink and an anti-kink in a chain. The photo- 
excited state is characterized by the relative momentum 
Qr = [Qe — Qh)/'^ between the momentum of the elec- 
tron pair Qe and that of the hole pair Qh, since the total 
momentum Qe + Qh is conserved in the photo-absorption 
processes. The several peaks in a{uj) shown in Fig. [5] are 
classified by qr- Based on this kink/anti-kink picture, we 
derive the effective model for the photo-induced dynam- 
ics in the V — t model in the limit oiV^t. We consider 
one kink/anti-kink pair. Details are presented in Ap- 
pendix |Al We obtain the transition probability for the 
current operator from the kink/anti-kink ground state 
with the relative momentum qri = tt/L to an excited 
state with qrf. The result is given by 



I {<lri,qrf) = I {qrf\j kri) \ =C — 



t I qrf sin qr 



It, 



(48) 



where a constant C — 144 for the lowest excited 
kink/anti-kink state with the momentum g^/ = qri + T^/ L 
and C = 64 for other g^/'s. The numerical results 
shown in Fig. [5] are well reproduced by this analytical 
method based on the kink/anti-kink picture. In partic- 
ular, it is clarified that the lowest energy peak around 
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FIG. 7: (color online) Optical absorption spectra in the ex- 
tended Hubbard model. Numerical data for rt =3, 6, 10, 15 
are almost overlapped. The inset is for the time dependence 
of superexchange- interaction energy defined by Eq. (|5ip . The 
parameters are chosen to be [/ = 20t, V = 5t, L = 9 and 
N = 5. 



Lo ^ 2. It marked by C corresponds to the transition 
given in Eq. (|48p with the minimum momentum transfer 
qrf — qri = 7r/i. lu the thermodynamic limit of i — > cxd, 
when the kink/anti-kink density n is fixed, a position of 



this peak is e 



and its weight is given by 



Dyt = fim nL ^ — 



9t 

= nL— ainqr 



(49) 



where £q is the energy of the spin-less fermion. Here, 
we neglect the interaction between kinks and anti-kinks. 
This corresponds to the Drude part of the optical con- 
ductivity. 

Some aspects in a{uj) in the DE model is explained 
from this charge kink/anti-kink picture. By changing the 
parameter values in the numerical calculation, we have 
confirmed that the main structures of a(w) in the DE 
model before pumping, e.g. the number of peaks, the V 
dependence of the peak position and others, correspond 
to the structures in the V —t model. Thus, we expect that 
the peaks in the DE model are classified by the relative 
momentum qr of the kink and anti-kink excitations, and 
that the photo-induced in-gap component corresponds to 
the Drude weight in the thermodynamic limit. On the 
other hand, the fine structures in aiuj) and the time- 
dependence of the in-gap component in the DE model 
are not explained by the kink/anti-kink picture and are 
sensitive to the parameter Jg. These are attributed to a 
time dependence of the spin structure as explained latter. 

As well as the V — t model, almost no-time dependence 
of the transient-optical absorption spectra is seen in the 
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extended Hubbard model in one-dimensional chain: 

+ rii-t-rijj. + nirij. (50) 

' (u> 

The numerical results of a{uj) are shown in Fig. [T] We 
adopt the one-dimensional finite-size cluster with the 
open boundary condition. The electron number per site 
is 0.5, and the pump photon energy and the dumping 
factor are set to be wq = 3.2t and 70 = 0.4i, respectively. 
As well as the DE model and the V — t model, the in-gap 
spectral component appears around uj ^ 2t hy photo- 
irradiation. This component does not show remark- 
able time dependence. We also calculate the exchange- 
interaction energy between the NN spins defined by 

^""' = Z7^E(^'-^^)' (51) 

with the spin operator Si (1/2) X^a/? 4a'''a/3Ci/3. It is 
shown in the inset of Fig. [7] that the change in E^^'° is 
less than 1% of the pump photon energy. Similar weak 
time-dependence in the spin sector was also observed in 
the half-filled Hubbard model<^ Through the analyses 
in the V ~ t model and Hubbard model, we conclude 
that the localized-spin degree of freedom which couples 
with the conduction electrons is essential for the time 
dependences of a{u)). 

C. Dynamical Correlations 

The transient electronic structure in the DE model is 
examined directly by calculating the one-particle excita- 
tion spectra. In Fig. [51 we plot A{q,uj) before and after 
photo-irradiation. The CO insulating gap is seen before 
pumping. The one-particle spectra are almost symmetric 
in the momentum space with respect to q = 7r/2. This 
is caused by the Brillouine-zone holding due to the stag- 
gered CO. There is no electron-hole symmetry in A{q, u) 
because of the localized spin degree of freedom. By photo 
irradiation, a photo-carrier band appears inside of the 
CO gap around w ^ 0, and spectral intensity of the up- 
per (lower) band around q = Q (tt) is weaken as shown 
in Fig.IUb) and[8|c). It is noticeable that a width of the 
photo-carrier band becomes broad with increasing time. 

The photo-induced charge and spin dynamics are ex- 
amined furthermore by calculating the dynamical charge 
and spin susceptibilities defined by Eq. and Eq. . 
Contour maps of the imaginary part of the dynamical 
charge susceptibility are plotted in Fig. [3] The exci- 
tation spectra before pumping show broad continuum 
around q = tt and sharp peak intensity at uj/t ^ 5.5 
and q ^ 7r/2. By photo- irradiation, excitation energy 
significantly decreases. Energy at the lower energy takes 




FIG. 8: (color online) The one-particle excitation spectra be- 
fore pumping and those at rt =3 and 10. Solid and dotted 
lines are for the electron and hole parts of A{q,iu), respec- 
tively. Momentum q of each panel increases from bottom to 
top. 

almost zero around q = and tt. Large spectral weight is 
seen around q = O.Stt, and total intensity is more than ten 
times larger than that before pumping. With increasing 
time, the band width of the continuum spectra increases 
weakly. 

Some features in N(q,io) are interpreted from the 
kink/anti-kink picture introduced previously. In the V — t 
model, which is equivalent to the = 1/2 XXZ model;^ 
excitation energy for a kink/ anti-kink pair is given by 

E{qt,qr) = V - Atcosqtcosqr, (52) 

in the limit of 1^ 3> t. The relative and total momenta are 
given by qr = {Qe-Qh)/'2 and qt = (Qe + Q/i)/2, respec- 
tively, where Qe(h) is the momentum for the NN electron 
(hole) pair. The total momentum corresponds to q in 
N{q,io). Before photo-irradiation, charge excitations for 
several relative momenta qr cause excitation continuum 
at fixed q{= qt) in N{q,Lu). By photo-irradiation, a pair 
of kink/anti-kink with momenta Qe and Qh{— ~Qe) is 
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FIG. 9: (color online) A contour map of the imaginary part of 
the dynamical charge susceptibility N{q, uj) before pumping 
and those at rt = 3 and 10. Data in the regions of g > 7r/2 
and q < n/2 are plotted in different scales. 



created. In particular, in the lowest photo-excited state 
obtained by tuning the pump energy at the optical edge, 
both Qe and Qh are almost zero. The excitation en- 
ergy in this case is E{0, 0) = V — At. In the excita- 
tion processes for N{q,ijj) in the one-photon absorbed 
state, ether Qe or Qh is changed, i.e. an excitation from 
(QcQh) — (0,0) to {q,0) or (0, g). Energy in this state 
is E{q, ±q) = V — 4t cos^ q. Thus, the charge excitation 



energy in the one-photon absorbed state is given by 

AE{q)^E{q,±q)-E{0,0) 

2<(1 - cos2q). (53) 

This momentum dependence explains the lower edge of 
the continuum spectra after photo-irradiation shown in 
Figs.Etb) and 9(c). 

On the other hand, continuum above the lower edge in 
N{q,uj) and increasing of the band width with increas- 
ing time are not explained only from this kink/anti-kink 
picture and are attributed to the localized spin degree of 
freedom. In the strong Hund-coupling limit of the DE 
model, the effective transfer integral teff between the 
NN sites depends on the spin structure in their sites. It 
is expected that by photo-irradiation, an average value 
of tef f increases with increasing time because of reduc- 
tion of the AFM correlation. This change in spin struc- 
ture contributes to the increasing of the band width in 
N{q,uj). 

Contour maps of the dynamical spin susceptibility are 
shown in Fig. [101 Negative intensity around w = is due 
to the fact that the initial state of S{q^uj) after pumping 
is not the ground state but the photo-excited state where 
the second term in Eq. (|25p is sometime larger than the 
first term. Before pumping, the spectral weight is seen in 
a region of the order oiuj ^ Js- We confimred a following 
dispersive feature in S(q,uj) by increasing Js (not shown 
in the figure); a sine-like dispersion in the lower energy 
edge, continuum spectra above the edge, and a strong 
intensity at g = tt. This feature is similar to that in the de 
Cloiseaux- Pearson mode in the one-dimensional spin-1/2 
AFM Heisenberg model. After the photo-irradiation, the 
characteristic dispersion disappears and the continuum 
spectra are broaden up to w ~ 2.5t. These changes imply 
a weakening of the AFM correlation by photo-irradiation 
and are consistent with the results of S{q) in Fig. |4l^b). 

D. Correlation between spin and charge dynamics 

In order to reveal the correlation between photo- 
induced charge and spin dynamics in more detail, we cal- 
culate the second moment of the in-gap band in A{q, ut) 
corresponding to the band width, the integrated spectral 
weight of the in-gap component in a(aj), and the correla- 
tion function between the NN localized spins. These are 
defined by 

W= f ^A(q,w)(a;-Wc)^<iw, (55) 
D= a{uj)dLo, (56) 
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FIG. 11: (color online) Time-dependences of the NN corre- 
lation between localized spins, the band width of the in-gap 
states in A{q, to), and the spectral weight inside of the optical 
gap in a{uj) denoted by Ks, W, and D, respectively. The 
data are subtracted by the data at rt — 3, and normalized 
by differences between their minimum and maximum values. 
We take lju ~ 5t, ul = —2.Qt, ui'jj = At and oj^ = Q.lt. 



pling between the charge and spin degrees of freedom in 
the photo-excited states. 

The time dependence of Ks a-t various values of a, 
which is a prefactor of the transfer integral in Eq. ([T]), 
is presented in Fig. [T^ a). The results are also plotted 
as a function of the time scaled by the transfer integral, 
rat, in Fig. [T2jb). The slope of the curve increases with 
decreasing a. By reploting as a function of time scaled by 
at, the slopes for several value of a are almost identical 
in a region of rt = 3 — 7. Through this analysis, we 
conclude that photo induced charge and spin dynamics 
are controlled by the electron transfer integral for the 
conduction electrons. 



FIG. 10: (color online) A contour map of the imaginary part 
of the dynamical spin susceptibility S{q,Lu) before pumping, 
and those at rt = 3 and 10. 



where uju and uil are the upper and lower edges of the 
in-gap component of A(q,uj), respectively, lo'jj and 
are those of the in-gap component of a{uj), and Wc is the 
center of mass for the in-gap band defined by 



(57) 



Numerical values are set to be uju = 5t, ujl = — 2.6t, 
co'jj = 4:t, and = O.lt. As seen in Fig. [Til all three 
curves show similar time dependence; intensity increases 
linearly after pumping and saturates around rt = 10. 
This identical time dependence indicates a strong cou- 



VI. PUMP-PHOTON DENSITY DEPENDENCE 

So far, we restrict our calculations to the one-photon 
absorbed states in a finite size system. In this section, we 
introduce pump-photon density dependence of the charge 
and spin dynamics in the extended DE model. Let us 
go back to the formulation for the electron-photon in- 
teraction. We consider a condition that amphtude of 
the pump photon is not weak and the interaction Hamil- 
tonian between electron and photon in Eq. ^ is not 
treated perturbatively. For convenience in numerical cal- 
culations, instead of Eq. ([7]), we introduce a Gaussian- 
type of the pump-photon vector potential given by 



^pump(T) = Ai exp 



21 



cosa;o(r - rg). 



(58) 
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FIG. 12: (color online) Time-dependence of the NN spin cor- 
relation between localized spins for several value of a. Hori- 
zontal axes are taken to be rt in (a), and rat in (b). 



where tq{> 0) is a center of the wave packet. We have 
checked that the results, such as the spin correlation 
function, calculated with this type of pump-photon vec- 
tor potential qualitatively reproduce the results with the 
Loretzian-type vector potential introduced in Eq. ([7|. 
Time evolution of the system is governed by a sum of 
the electronic part and the interaction part, i.e. 'H(r) = 
"Ho -I- T-L'ij) where Ho and 'H'(r) are given in Eq. ([T]) 
and Eq. respectively, with Eq. ([55)) . We assume that 
To ^ 7(^^, and the system at time t = is in the ground 
state of Ho, i.e. |(/>(t = 0)) = |0). Time evolution of the 
wave function is calculated by Eq. p9| from \(j){T = 0)). 
In the Lanczos method, we confirmed the parameters 
M2 = 20 and 5t — 0.01/t are sufficient to obtain the 
reliable results. In the multi-target DMRG algorithm, 
we adopt {|0(t)) , |(/)(t + dr))} as the target states in the 
mixed-state density matrix defined in Eq. pop . We chose 
the parameters M2 = 20 and 5t = 0.1/t in Eq. (pg)) . and 
the truncation number m = 600 — 800. In Eq. ([55)l . the 
pump-photon energy is chosen to be ujq = bt which is 
around a middle of the lowest absorption band shown in 
Fig. [3] We chose the damping factor 70 = 2t in order 
to demonstrate clearly the pump-photon density depen- 
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FIG. 13: (color online) Time dependence of the absorbed pho- 
ton density nph(j) for several value of (solid lines), and 
that of the pump-photon vector potential Apump(T) in the 
case oi A\ — 1 (broken line). 



dence. System size is chosen to be L =9, 13 and 17. 

In Fig. fTnT a). the photon densities absorbed in the sys- 
tem are plotted for several values of a prefactor A\ in the 
pump-photon vector potential in Eq. (|58p . We define the 
absorbed photon density at time t by 



nphir) 



(H(r))-(H(-^)) ^ 
LujQ 



(59) 



where H{—oo) is the Hamiltonian at time t <C tq. We 
also plot ^pump(''') in the same figure. The absorbed 
photon density increases abruptly around r = tq = 3/i, 
and is saturated around t = 5/t. A saturated value of 
Uph monotonically increases with increasing A^. 

We focus on the electronic states at r = 5/< where 
^pump(''') is almost fully damped and nph(T) is saturated. 
The following results are not sensitive to a choice of a 
value of r. In Fig. [Tif a). the absorbed photon density 
and the charge and spin correlation functions at time 
T = 5/t are plotted as functions of the pump-photon 
intensity. We define that AKjyf^s-j is a difference of Ks{N) 
at T ~ 5/t from that at r <C tq. The spin correlation 
function Ks was introduced in Eq. ([M)) and the charge 
correlation function is defined by 



Kn = 



{•y> 



(60) 



Both the two increase with increasing A^ and tend to 
be saturated around = 3 — 4. In order to compare 
the spin and charge dynamics in more detail, we plot the 
first derivatives of AKs and AKn with respect to A^ in 
Fig. [Hfb). It is shown that a change in AKs is more 
sensitive than that in AK^ in weak photon density, and 
is almost saturated around A^ = 3. 

To clarify this different spin and charge dynamics fur- 
thermore, we plot changes in the spin and charge corre- 
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FIG. 14: (color online) (a) The absorbed photon density 
nph{T), the charge correlation function (L — 1)AKm, and 
the spin correlation function (L — l)AA's as functions of 
the pump-photon amplitude, (b) The first derivative of Km 
and Ks with respect to Al, i.e. (L — l)d{AKN) /d(Al) and 
(L - l)d{AKs)/d{Af). Time is chosen to be tr = 5. 

lation functions divided by the absorbed photon number 
Nph as functions of Uph in Fig. [151 We introduce 

Nph = Luph. (61) 

Values of AK]\j /Nph are almost independent of Uph and 
are identical to be about one. This result does not 
show remarkable size dependence. On the other hand, 
AKs/Nph shows clear nph dependence and almost lin- 
early decreases with increasing Uph- A weak system size 
dependence is seen in AKs/Nph in a region of small Uph- 
This difference in the spin and charge dynamics is in- 
terpreted as follows. The result of AK^/Nph ~ 1 im- 
plies that one-photon generates one electron-hole pair. 
This is explained by the kink/ ant i- kink picture in a one- 
dimensional chain introduced in the previous section. On 
the contrary, it is shown that destruction of the AFM 
spin correlation by pumping is effective in the weak pho- 
ton density case. This may be attributed to the inter- 
action between the photo-carriers. In the case of large 
V , the photo-carriers are not exchanged with each other 
in a chain. Each carrier destroys the AFM spin corre- 
lation through the DE interaction. In the case of high 
photo-carrier density, i.e. large 7ip/i, regions where photo- 
carrier destroys the AFM correlation are overlapped and 
destructions of the AFM correlation are interrupted with 
each other. 



^(L-1)AKN/Nph 
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FIG. 15: (color online) Differences in the charge and spin 
correlation functions devided by absorbed photon number. 
Data are plotted as functions of the absorbed photon density. 
Solid, broken and dotted lines are for the results in system 
sizes L =9, 13 and 17, respectively. 
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FIG. 16: (color online) The kinetic energy as functions of the 
pump-photon intensity. Broken line is for the results in the 
spin-less V — t model. 



Roles of spin degree of freedom in photo-induced dy- 
namics are also seen in the photon-density dependence 
of the kinetic energy. In Fig. 1161 the kinetic energies in 
the extended DE model and the V ~t model are plotted 
as functions of the pump-photon intensity. In the cal- 
culation of the V — t model; one dimensional cluster of 
L = 9 and N — 5 with the open-boundary condition is 
used. The parameter values are chosen to be 1^ = 5t, 
LjQ = 5t and 70 = 2t. We plot Et defined in Eq. ([39]) and 
its correspondence in the V ^ t model at t — 5/t where 
Uph is saturated. A non-monotonic dependence on is 
seen in the DE model; Et decreases at first in a region 
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of weak photon density, and increases above ~ 0.5. 
This is in contrast to the results in the V — t model where 
Et monotonically increases with the photon amplitude. 
That is, the spin degree of freedom qualitatively changes 
the photo-induced dynamics in a weak intensity region. 
This observation is explained from a competition between 
carrier excitation by photo irradiation and the coopera- 
tive weakening of the charge correlation and the AFM 
correlation. A monotonic increasing of the kinetic en- 
ergy is caused by a carrier excitation by pump photon. 
On the contrary, the charge correlation is reduced by a 
weakening of the AFM correlation due to photo-carrier 
motion. In a region of small Af, the latter overcomes 
the former, and, and as a result, Et decreases with Af. 
On the other side, in the strong pump-photon density, 
a destruction of AFM correlation is saturated, and the 
photo-carrier excitation effect emerges. 



VII. DISCUSSION AND CONCLUDING 
REMARKS 

In this section, we summarize the present numerical 
results and discuss implications for the photo-induced 
phenomena in perovskite manganites. The present main 
results are listed: 1) a new state appears inside of the 
insulating gap by photo- irradiation (see Figs. |3] and 
2) time dependence of the NN spin correlation Ks is 
strongly correlated with those of the band width W of 
the in-gap states in A(q, co) and the spectral weight D 
inside of the optical gap in a(w). These are scaled by a 
universal curve (see Fig. [TT|) . 3) The characteristic time 
scale Ts, when a change in Ks is saturated, is governed 
by the electron transfer t (see Fig. [12]), i.e. ts ^ 10 /t. Al- 
though the present our calculations have been carried out 
in one-dimensional finite-size clusters, the characteristics 
in the photo-induced dynamics listed above arc quali- 
tatively reproduced by the time-dependent unrestricted 
Hartree-Fock method in two-dimensional systemsi^ By 
using the realistic parameter values for manganites, this 
time scale ts is about few fs. The present calculations 
predict that both the CO correlation and the short-range 
AFM correlations collapse cooperatively within this time 
scale and a metallic state appears in the optical spec- 
tra and the photoemission spectra. This ultra-fast CO 
melting and generation of a metallic state are consistent 
with the recent pump-probe experiments in manganite; a 
melting of charge order within about 30fs by irradiation 
was observedji^ 

We note that the weakening of the NN spin correlation 
within few fs is not contradict to the magneto-optical 
Kerr spectroscopy measurement. In several charge or- 
dered manganites, it was shown that the photo-induced 
spin dynamics observed by the Kerr spectroscopy is 
slower than the photo-induced charge dynamics probed 
by the optical reflection i^^i^^i^^i^^i^^ These experiments 
and the present theory suggest that the short-range and 
long-range spin dynamics are governed by the different 



mechanisms, the electron transfer and the spin-orbit cou- 
pling, respectively. The macroscopic spin dynamics can- 
not be dealt with by the present ED method where the 
spin angular moment is conserved in time evolution. This 
issue was studied in detail in our previous paper— by 
utilizing the time-dependent Hartrec-Fock method with 
spin relaxation processes. Photo-induced dynamics in 
the macroscopic magnetic moment is controlled by the 
phenomenological spin relaxation constant F, and the fi- 
nite magnetic moment appears in the longer time scales 
of Ti ~ 1/r than Ts- One noticeable point is that change 
in the charge sector from the CO to the metal is almost 
completed within the short time scale ts , and growing of 
the macroscopic magnetic moment around tl docs not 
associated with change in both the short- and long-range 
charge correlations. This is attributed to the fact that the 
spin-orbit interaction, by which the conservation of the 
macroscopic spin angular moment is only broken, play a 
crucial role on the photo-induced macroscopic spin dy- 
namics. 

It is shown that, even in the short-time scale, photo- 
induced spin and charge dynamics does not always show 
identical behaviors; the change in the short-range spin 
sector becomes remarkable in the case of weak pump- 
photon density, in comparison with that in the charge 
sector. This is seen in the spin correlation function and 
the kinetic energy as functions of the pump-photon den- 
sity (see Figs. 1141 and [TB)) . These different dynamics be- 
tween charge and spin are caused by the fact that the ini- 
tial CO order is collapsed directly by the photon, while 
the AF spin correlation is collapsed by carrier motion. 
The detailed measurements of the total weight of the op- 
tical conductivity spectra may detect the non-monotonic 
behavior of the kinetic energy as a function of pump- 
photon density. 

In conclusion, we study numerically the photo-induced 
electronic dynamics in the generalized DE model in one- 
dimensional chains. The Lanczos and DMRG methods 
are utilized to calculate the time evolution after photo- 
irradiation and the several spectral functions. Roles of 
the localizcd-spin degrees of freedom in the DE model are 
focused on. By pump-photon irradiation into the CO and 
AFM insulating state, the collapsing of the AFM corre- 
lation and the appearance of a metallic state occur coop- 
eratively. This time evolution is governed by the electron 
transfer integral of the conduction electron. The numer- 
ical results are explained by the charge kink/anti-kink 
picture and their spin-dependent kinetic motion. The 
pump-photon density dependence of the spin and charge 
dynamics are also examined. Roles of the spin degree of 
freedom are remarkable in the case of the weak pump- 
photon density. 



Appendix A: V-t model in one dimension 

In this APPENDIX, the effective model obtained from 
the V — t model and an analytical formulation for the op- 
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tical absorption spectra in this model are presented. We 
start from the one-dimensional spin-less V — t model in- 
troduced in Eq. (|47|) where the numbers of sites and elec- 
trons are L and i/2, respectively. In the limit of F ^ i, 
the groimd state electron configurations are denoted as 

• • ■ 101010101 • • • where 1 and represent electron occu- 
pied and unoccupied sites, respectively. The low energy 
excited states arc given by the configurations where one 
NN electron pair and one NN hole pair exist shown as 

• • • 100101101 • • • . In this subspace, the effective Hamil- 
tonian in the limit of ^ t is obtained by the canonical 
perturbational expansion. In the case where the electron 
pair locates in the left side of the hole pair in a chain, 
the Hamiltonian is given as 



L/2 



vt 



f ajflj + b\bi+i + h.c 



We introduce the two fermion operators, 

bi = {1- n2i)4i-i - (1 - n2i-i)dl^ 



(Al) 



(A2) 



(A3) 



which correspond to the annihilation operator of the NN 
electron pair and that of the NN hole pair, respectively. 
We impose the relations 



(A4) 



and wc have the following anti-commutation relations 

{aj, aj = Sij {n2i + n2i+i) , (A5) 
{blbj} = % [(1 - n2^) + (1 - n2.+i)] , (A6) 

and other anti-comutators are zero. The operators in 
Eq. ((J2|l and Eq. satisfy the conditions 



0. 



(A7) 



In this model with the open-boundary condition, the 
eigen states are classified by two momenta ki and k2 for 
a electron pair and a hole pair as 



|fci,fc2) =— — - ( sin Xjki sin Xj fc2 ~ smxik2smxjki 

X ajfe] |0) . (A8) 
The corresponding eigen value is 

Ek,,k2 = -2t {cos ki + cos k2) . (A9) 



We introduce the momentum k = mr/{L + 1) and an 
integer number n(= 1, 2, • • • , i). In the present scheme, 
the current operator is given by 

J = -it E - h\b,+i - H.c}j . (AID) 

i 

A matrix element of j between two states is obtained by 



2 

it I dxj 



dxi 



{kfikf2\j \ kiih2) 

\^/ Jo Jo 
X [s'mxikfi smxjkf2 — s'mxikf2 sinxjkfi) 

X 2sinfcii {cosxikiisinxjki2 + cos Xiki2 sinxjkn) 

— 2 smki2 (sinxikii cosxjki2 -\- su\Xiki2 cos Xjkn) 

(All) 

In the thermodynamic limit L — s- 00, the total momen- 
tum fc*°'(= fci — ^2) is zero in the optical process, and 
the eigen states are characterized by the relative momen- 
tum q[= {ki + ^2)72]. By carrying out the integrals in 
Eq. (IA11[) . the transition probability from the state with 
Qi to that with g/, i.e. I{qf,qi) = is given by 



62^2 



(A12) 



for the minimum momentum change qf = + tt/L, and 



Hl^^Qf)-- ^2 



r / g/sinq. 
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(A13) 



for others. Above two are summarized in Eq. (|48 
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